DTP 95-75 
April 1996 



o 



Topological Chern-Simons vortices in the 

0(3) cr-model 



on 

J. GladikowskiJ 

Department of Mathematical Sciences 
. South Road, Durham DH1 3LE, England 

OO 

^SJ ' Abstract 
> 

We present a classical, gauged 0(3) cr-model with an abelian Chern-Simons term. It 



shows topologically stable, anyonic vortices as solutions. The fields are studied in the case of 



rotational symmetry and analytic approximations are found for their asymptotic behaviour. 
The static Euler-Lagrange equations are solved numerically, where particular attention is 
paid to the dependence of the vortex' properties on the coupling to the gauge field. We 
compute the vortex mass and charge as a function of this coupling and obtain bound states 
for two-vortices as well as two-vortices with masses above the stability threshold. 

> : 

1 Introduction 

a; 

The 0(3) cr-model in 2-dimensional Euclidean space is a classical field theory which sup- 
ports soliton solutions [Q] . Its scale invariance can be broken by the addition of a potential 
term. This does not prevent the soliton from shrinking, however, its size can be fixed by 
the inclusion of higher order terms in the field gradient Q . An example for such a theory 
is the baby Skyrme model ||. Alternatively, the scale invariance of the 0(3) cr-model can 
be removed and the soliton be stabilised (at least in principle) by gauging a (7(1) subgroup 
of the fields internal symmetry group ||, §|. The dynamics of the U(l) gauge field in 
such models is ruled by Maxwell- and/or Chern-Simons actions. For each of these cases 
potential terms have been constructed such that the corresponding models yield self-dual 
equations of Bogomol'nyi type. The potential term determines the asymptotic behaviour 
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of the fields which can either obey the gauge symmetry & 0] or break it ||. The models 
with broken gauge symmetry show topological soliton solutions with quantised magnetic 
flux and in this they resemble the well-studied vortices in the abelian Higgs model and its 
generalisations , see for instance ||, . 

Here we investigate a gauged 0(3)cr-model with Chern-Simons action. Chern-Simons 
theories are an object of intense research because their quantised version is relevant for sys- 
tems of strongly correlated electrons e.g. in superconductors or in the quantum Hall effect 



[11]. In this paper we consider a static classical Chern-Simons model, whose potential term 
preserves the gauge symmetry and is chosen to produce exponentially localised configura- 
tions. They carry fractional angular momentum and have a lower topological bound on 
the energy which is, however, not saturated. We solve the equations of motion numerically 
for radially symmetric fields and study the dependence of the solutions on the coupling 
strength to the gauge field. We also look at two vortices on top of each other and on their 
mutual attraction dependent on their coupling. The asymptotic behaviour of the fields is 
studied analytically and conclusions about intervortex forces are drawn. 

Recently, static solitons were found in a gauged CV 1 model which includes a Chern- 
Simons term and a potential term equivalent to the one considered here |l2| . In its standard 
version the CV 1 model represents merely a different choice of fields to the 0(3) cr-model. In 



[12], however, the gauged symmetry is the internal U(l) symmetry of the two-component 
complex CV 1 vector which lies on S 3 . Therefore we expect our solutions to be different to 



the ones presented in [12], but it is nevertheless instructive to compare them. 



2 Chern-Simons solitons revisited 

We consider the following Lagrangian of a gauged 0(3) cr-model in (2+1) dimensions. It 
contains a potential term and the behaviour of the gauge field A a is governed by a Chern- 
Simons term 

L = \ {D a 4>f - ^d a A^ - v\l - n • 0) . (1) 

The fields </> are three-component real vectors and subject to the constraint (f> - (f> = 1, hence 
they take values on the two-sphere S^. The metric is chosen to be g a p = diag(+, — ,— ). 
Throughout this paper, greek indices run from to 2 while Latin indices denote the two 
spatial dimensions 1,2. We work in geometrical units in which the velocity of light c = 1. k 
and /i are real coefficients of dimension length and 1/length respectively and for dimensional 
reasons the Lagrange density (|l]) should be thought of being multiplied by an overall factor 
of dimension energy. We fix our mass scale by putting this factor to one. The fields (f> will 
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be frequently referred to as matter-fields (in distinction to the gauge fields) and to S| as 
the iso-space. The potential term in ([!]) reduces the symmetry of the model to O(2) iso , i.e. 
to rotations and reflections perpendicular to the vector n. It is this symmetry that is to 
be gauged and by choosing n = (0,0, 1) we select the SO(2) iso subgroup which consists of 
unimodular rotations about the z-axis. D a (<f>) is the covariant derivative and given by: 

D a cj) = d a cj) + A a (n x cj)) . (2) 

The ungauged Lagrangian shows symmetry under combined reflections in space and iso- 
space: 

P : (x 1 ,x 2 ) ^ (-xi,x 2 ) and C : (0 X , fa) ->• (-01, fa) > ( 3 ) 

which can be thought of as a parity operation and charge conjugation. The Chern-Simons 
term breaks the parity symmetry explicitly by changing its sign under P. It also breaks the 
time-reflection symmetry T which corresponds to Aq — * —Aq. However, the Lagrangian is 
still symmetric under CPT. 

The potential term can be thought of physically as an analogue to the Zeeman coupling 
between spin fields and an external, constant magnetic field in n-direction with coupling 
strength /j, 2 . Such terms occur for example naturally in the description of the quantum Hall 
effect. 

Because we are interested in configurations with finite energy, we require that the po- 
tential term and the covariant derivative vanish at spatial infinity. Hence we impose: 



lim d>(r) = n . (4) 

r— »oo 

This boundary condition allows to one-point compactify the physical space such that fields 
are maps: 

: Sl^S 2 . (5) 

These maps are elements of homotopy classes which form a group isomorphic to the group 
of integers. This integer or degree N counts the number of times 5"? is covered by a single 
covering of S 2 . It can be written as the integral over the zero component of the topologically 
conserved current: 

la = 77-£a/3 7 " (9/30 X <9 7 0) , (6) 

such that the degree N is obtained from 
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N= / d 2 xl Q , (7) 



where the range of integration is S^. Following the standard convention we will call finite 
energy solutions with N > vortices and those with N < antivortices. (^) is in fact the 
topological charge of the 0(3) a fields and it is not obvious that this is also a topologically 
conserved quantity in the gauged model. We will therefore address this question again 
below. 

The equations of motion derived from (||) can be written in terms of the matter-current 
J a and the electromagnetic current j a 

J a = 4> x D a 4> , j a = n-J a . (8) 

The Euler-Lagrange equations are 

D a 3 a = /i 2 (n x 4>) (9) 
d p A^ . (10) 



Note that by (|l0|) the gauge fields are completely determined by first order equations and 
do not have own dynamics in the strict sense. Equation ([!(]) for a = is Gauss' law 

kB(tcy x (b) , . 

where we have used that n = (0,0, 1) and B = eoijd l A^ , taking £012 = 1- The equation of 
motion ( |io[ ) implies that for non-singular A a the electromagnetic current j a is conserved 
(d a j a = 0)- The current can be written conveniently as j a = (p,ji), where p is the charge 
density of the soliton while ji denotes its electric current. The Lagrangian (jlj) can be 
expressed in terms of j a : 

£ = 1 {D a <t>f - X -A a f + p\\ -n-4>). (12) 

This shows explicitly that the gauge fields A a are coupled to the electromagnetic current 
j a . The electric field E and the magnetic field B are related to j a as follows: 

P TP . J 1 



5 = -- , E i = e ij J -. (13) 



K 



The first equation leads to a relation between the magnetic flux $ and the electric charge 
Q of the configuration 
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$= / d 2 xB = -- [ d 2 xp = --. (14) 

J K J K 

The theories energy-momentum tensor is obtained by the variation of the Lagrangian with 
respect to the metric g a p 

T aP = {D a cf>){D^) - g a(3 Q(D 7 0)(Z^) - n\l - n • 0)) . (15) 

The integral over the component Tqq is the total energy of the soliton 

E cs [<j>,A] = |d 2 xi( J Do0) 2 + ^(A0) 2 + / u 2 (l-n-0). (16) 

Note that the Chern-Simons term does not contribute directly to the energy because of 
its metric independence. The rotational symmetry of the Lagrangian leads to a conserved 
angular momentum M of the soliton 

M = J d 2 x (x x p), (17) 

where the cross product stands for X1P2 — P\X2- M is a vector pointing perpendicular out 
of the plane of motion. The components of the momentum density p are given by 

Pi = T 0l = D <f> • Di<f> . (18) 



3 A bound on the energy 

Next we give a proof that E cs [cj>, A], the energy in our model given by (|l~6|), is bounded from 
below by a topologically conserved quantity. This is not obvious, because the gauged pure 
0(3) cr-model does not have a lower bound on the energy, unlike its ungauged counterpart, 
where the solutions saturate the Bogomol'nyi limit. The first step in the proof is to use an 
auxiliary energy functional i? aux [0, A] which is of Bogomoln'nyi type and was constructed 
in ||]. Because we wish this section to be self-contained, we will repeat below parts of the 
analysis given in this reference. First, we show that the energy gap between E cs and S aux 
(or a multiple of it) is positive and then complete the argument by demonstrating that 
£ mi > 47r|JV|. £ aux reads as: 

£ aux [0, M = \J d 2 x B 2 + {D^f + (1 - n • </>) 2 . (19) 
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In order to be consistent in the notation of dimensions, both the potential term and the 
magnetic field must be thought of being multiplied by a parameter of dimension 1/length 
squared and length respectively. These parameters are of magnitude one and subpressed 
in (|l9|). To compare E cs with _E aux one first observes that (Dq<P) 2 > k 2 B 2 , due to Gauss' 
law (pA|). Now we carry out a rescaling of x in our functional E cs , namely x — > kx, which 
transforms B — > B / k 2 and 4>(x) — > (f)(nx). The potential term then reads as k 2 /j, 2 (1 — n • (p) 
and is greater than (1 — n • (f>) if k > l//i. To verify that E aWK is smaller than E cs we use 
that since < (1 — n • <fi) < 2, it follows that (1 — n • 0) > 4(1 — n • (f>) 2 and one sees that 
for £7 aux holds 

£cs>£a UX if «>l//i. (20) 

In the case k <1//j, we assess an energy bound by multiplication of each individual term in 
the energy density with n 2 fj?. This gives 

E cs > kV#«« if K < 1/n . (21) 

This already proves the bound for 2? cs , but it is instructive to see in detail that I? aux defines 
a Bogomol'nyi model. In order to achieve this, we rewrite the auxiliary energy functional 
as 

E^[B,cj>] = ^j d 2 x (D 1( f>±cf>x D 2( f>) 2 + (BT(l-n^)f ± J d 2 xL . (22) 

Lo is composite of the cross terms and can be understood as the zero component of the 
solitons gauge invariant topologically conserved current: 

L a = e aPl (0 • {Df>4> x + dPAr(l - n • 0)) . (23) 

Up to a surface term, this current is equivalent to l a , the topological current of the ungauged 
0(3) cr-model (|6[). If the solutions are required to have finite energy, then must tend to 
zero faster than 1/r as r goes to infinity, hence it follows by Stokes' theorem that the surface 
term integrates to zero. It was pointed out in || that the conserved topological charge 
(integral over Lo) equals the degree N of the map cj) if the gauge symmetry is unbroken 
(as it is in our case) but can differ from N for broken symmetry. In order to saturate the 
Bogomol'nyi bound, both squares in (22) have to vanish, such that the following two (anti-) 
self-dual equations are read off 

Di0 = T0x£»20 , £ = ±(l-n-0). (24) 
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These equations were discussed in detail for a special choice of the fields in ||. There it was 
shown that they yield a one-parameter family of solutions which are degenerated in their 
energy but differ in their magnetic flux. 

By using the sign ambiguity in front of the integral over Lq in (|2^) we can restrict our 
discussion to the case B > and the upper sign without a loss of generality. Equation ( p2[ ) 
then implies 

£ aux > / d 2 xL = 4ir\N\ , (25) 
The equality holds for self-dual solutions. 

4 Static vortex solutions 

To find static solutions in our model we restrict ourselfs to the two-dimensional hedgehog 
[13], which is in terms of the polar coordinates (r,9): 



( 



4>(r,0) 



sin/(r) cos nO ^ 
sin/(r) smn9 
\ cos/(r) j 



(26) 



For this field the topological charge density, the integrand of (Q), equals 

71 i 

Z ° = 4^ /Sin/ ' (27) 

where the prime denotes the derivative with respect to r. By integration one easily sees 
that n = —N. For the gauge field A a , the most general ansatz which leads to radially 
symmetric and static fields is given by 

A = nv(r) , A e = na(r) , A r = h(r)t , (28) 

where t denotes the time and the factor n is introduced for convenience. We fix our gauge 
by putting A r = and obtain the equations 



, / 2 /> + ir 2 



/ +— = n 1 i t J -- v sin/cos/ + /i 2 sin/ (29) 

r \ r z I 

1 (a + 1) . 2 

v = sin j . (30) 

K r 



Gauss' law ( pTT] ) reads in terms of a, v and /: 
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a = rv sin 2 / . (31) 

K 

We are interested in finite energy configurations, which requires that D a {4>) — > as r — > oo. 
To guarantee this and the regularity of the fields at the origin we impose the following 
boundary conditions 

o(0) = 0, /(0)=7T, v(0)=v 

lim^oo f(r) = , linir^oo a(r) = , lmv^ v(r) = t) m 

where vq, Vqq and are constants. With these boundary conditions it is clear that constant 
fields a and v are not a solution of ( |30j ) and (|3l|), which can be shown by contradiction. If 
a were a constant it would have to be zero everywhere because of (|32|) , in which case ( |30|) 
implies that i> is not a constant which in turn, via (^Tj) leads to a non-constant a. A similar 
argument applies for the case of v being constant. Hence the Euler-Lagrange equations do 
not lead to vanishing flux and charge. 

The total (static) energy is given as the integral over the energy density e, which reads 
in terms of the fields /, a and v (Jltf ) 



f n 2 ({a+lf 



+ v 2 sin 2 / + ,u 2 (l - cos /) . (33) 



2 2 \ r 2 j 
For the angular momentum (|l7]) one obtains 

M = -TTKNa^aoo + 2N)n . (34) 

Hence one sees that the angular momentum of the vortex is fractional and the vortices are 



(classical) anyons. The electromagnetic fields (13) are radially symmetric by construction 
and read as 

B = N— , E r = Nv . (35) 
r 

The electric charge and magnetic flux are not topologically quantised (unlike in the abelian 
Higgs model, for instance) and depend on the parameters in the model 

$ = N [ rdrdO — = -2-nN aoo = -— . (36) 
J r k 
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5 Asymptotics 



The boundary conditions (32) allow us to derive asymptotic approximations to the equations 
of motion (^). By approximating sin / ~ / and cos / ~ 1 for large r, the equation for / 
simplifies to 

f 1 ((„ _L^2 \ 

(37) 




where 

k 2 = fj 2 - r^vlc . (38) 

The asymptotic solution of / depends on the value k takes. There are three possible cases. 

1. ) \fi\ > | m>oo|) k real. 

The solution to ( j37|) for real k are given by modified Bessel functions / ~ K m (kr),m = 
^(floo + 1) with the asymptotic behaviour 

f~^e~ kr , (39) 

This shows that k can be understood as the effective mass of the matter fields (f>, by 
denoting the inverse decay length of the profile function. The asymptotics of the field are 
determined by the potential term which defines the theories vacuum structure. Therefore it 
is not a surprise that the vortex' matter field looks asymptotically like the baby Skyrmion 
investigated in Q, where the same potential term was used. 

2. ) \i < |n«oo|, k imaginary. 

This case leads to oscillating fields with an amplitude that falls off proportionally to 1/ yjr 
in leading order. The substitution k = ik in fl39| ) verifies this instantly and also shows that 
k is proportional to the inverse wavelength of the oscillations. The energy density of these 
solutions behaves asymptotically like \jr in leading order and hence the energy of these 
fields is infinite. This is, of course, not a physically relevant solution such that we exclude 
it from our further discussion. 

3. ) /i = | m;^, k = 0. 

The critical case is in fact just a special case of 1.), with vanishing exponential such that 
the profile function / ~ 1/r. The energy of these solutions is also infinite, because the 
leading term in the energy density is proportional to v^f 2 . Numerically we find that all 
these solutions occur but restrict our discussion to the case 1.), which gives the following 
constraint on the solutions: 

| A* | > \nVoo\. (40) 
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Using expression (^) , we find for the electric and magnetic field in the limit of large r 

t-t 1 —2kr d 1 —2kr f A-i\ 

hi r ~ — e , is ~ — e . (41) 
r r 

This shows that the electromagnetic fields fall off much faster than the matter field /. 
Therefore the electromagnetic interactions are expected to be negligible in the context of 
long-range vortex interactions. The electric field is a vector lying in the plane of motion 
while the magnetic field B can be thought of as pointing perpendicular out of the plane of 
motion. Its asymptotic shape is similar to the one of the Skyrme-Maxwell soliton discussed 
in ||, where it was argued that such a magnetic field resembles a magnetic dipole in two- 
dimensional electrodynamics. 

For small r the fields can be approximated by power series 

f^ir + cr^, v^vo + dr 2 ^, a « gr 2 ^ +2 , (42) 

where c and vq are free parameters while d and g are given as functions of n, k, c and vq. 
Note that for finite energy solutions c and vq are not completely independent on each other. 

For the Skyrme-Maxwell solitons it was found that the electromagnetic short range 
interaction decreases the energy per soliton and in particular leads to more strongly bound 
two-soliton states. Here, having a non-zero electrical charge distribution we expect this 
effect to be weakened by the Coulomb repulsion of the solitons electric field. 

6 Numerical results 

We solved the set of equations ( p9|) numerically by using a shooting method and a relaxation 
method. For both the shooting method and the time evolution in the relaxation method 
we employed a fourth-order Runge-Kutta method. In order to perform the numerical 
integration we had to fix the parameters in our model. Using geometric units in which 
the energy and length are of unit one, we are left with [i and k to be fixed. However, 
the parameter space is in fact one-dimensional which can be verified by carrying out the 
rescaling x — > kx, B -> B/k 2 . Thus we can fix [i for all our computations without a loss of 
generality. We choose fx = vO.l, a value which allows us to compare our numerical results 
with the ones obtained in the Skyrme-Maxwell model || where the same value has been 
used. 

We looked at the dependence of the solutions of degree N = 1 and N = 2 for a range of 
K. This parameter determines the strength of the Chern-Simons term and is proportional 
to the square root of the inverse coupling to the gauge field, which can be seen by a simple 
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substitution A a — > A a /i/K. 

Fig. la) shows the dependence of the static energy or mass on k. Small k which corresponds 
to strong coupling leads to lighter vortices for both the one-vortex and the two-vortex. For 
large k the mass or static energy E cs tends to a constant but remains relatively close to 
the Bogomol'nyi bound, staying below 1.1 (in units of 47r|iV|) for the one-vortex and the 
two-vortex. Thus our vortices are significantly lighter than the gauged baby Skyrmions, 
which tend to a mass of E SM = 1.546 for weak coupling. The energy gap arises partly due 
to the Skyrme term which is not present here. 

It is particular interesting to look at the relative static energy per vortex. We denote 
the energy of the one (two) vortex by E 1 (E 2 ) . The energy difference AE = E 2 — 2E 1 can 
be interpreted as binding and excess energy of the two-vortex for AE < and AE > 
respectively. In the case AE < the vortices form bound states while for AE > we 
expect that vortices on top of each other are unstable under perturbations and experience a 
repulsive force. From Fig. la) it is clear that in our model both cases occur. For small k the 
two-vortex is in an attractive regime as it is for large k, however, there is an intermediate 
region n l ci < k < for which the two-vortex is unstable (in the sense that its decay is 
energetically favourable). k ci is the critical coupling for which AE = 0. Numerically we 
find that k 1 ci = 0.632 and < r = 2.215. 

This result can be explained in a semiqualitative way. In the limit of large k the gauge 
fields decouple from the matter fields and become very small as compared to the matter 
fields. The study of ungauged solitons (e.g. in Q) showed that pure matter forces are 
often attractive for two-solitons. This is also the case in our model. For very small k, 
however, the magnetic flux tends to a constant (the reader finds the explanation for this 
below) while the repulsive Coulomb force is oc Q 2 oc <5 2 ac 2 , thus becoming weaker with 
increased coupling, see Fig. lb). We found numerically that increased coupling leads to a 
stronger bound configuration and in this our model is similar to the Skyrme-Maxwell model. 
The intermediate or repulsive range can be understood as a regime in which the Coulomb 
repulsion dominates the attractive forces of matter and magnetic field. It is within this 
range that the electric charge has its maximum value Q max = Q(K max ), where numerically 
«m« = 0.75 (for N = 1) and K max = 0.92 (N=2). 

In showing such a behaviour, the vortices resemble the fields of the abelian Higgs model 
where a similar transition between repulsive and attractive regime occurs, depending on the 
strength of the potential term. This characteristic is used to describe temperature driven 
phase transitions between type-I and type-II superconductors. The shape of the vortex 
shows a strong dependence on k, which is foreseeable by the interpretation of n as the 
coupling parameter to the gauge field. Because the vortex' Coulomb interaction is repulsive 
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it favours a spreading of the soliton. In agreement with this picture, we obtain vortices with 
their maximum width range at k ~ 1, where the electric charge takes its maximum values. 
See also Fig. 2a) and 2b), where the energy density e and the profile function / are plotted 
in dependence on the coupling. 

On the other hand, if the electromagnetic interaction is coupled only very weakly , 
the Lagrangian reduces to the 0(3)<r-model term plus the potential term, such defining a 
configuration which is known to be unstable against shrinkage due to the Hobart-Derrick 



theorem [14|. In accordance with this discussion one sees from Fig. 2 that for large k the 
vortex becomes more localised. This clearly shows that the potential term in the Lagrangian 
favours a shrinkage of the vortex. 

We also observed that increased magnetic flux as it occurs for small n leads to more 
localised solitons, like it does for Skyrme-Maxwell solitons Q. For very small k both the 
gauge fields a and v tend to singular configurations at the origin. In this limit, the gauge 
field v takes large values at the origin and is zero everywhere else, while a tends to — 1 
everywhere except at r = 0, which is fixed by its boundary condition, see Fig. 2c) and 2d). 
In that the behaviour of a is similar to the gauge field in the Skyrme-Maxwell model. We 
conjecture that the origin of this coincidence is the particular ansatz chosen for the gauge 
field, which leads to terms in the energy density depending on (a+1) 2 and so makes the value 
Oqo = — 1 exceptional. The strong coupling limit therefore leads to dynamically quantised 



flux and in addition implies via (36) that the electric charge vanishes for k — > 0. For the 
electric and the magnetic field we find that they form a ring (cf. Fig. 3b), a feature which 
was also observed by Jackiw and Weinberg in a self-dual Chern-Simons model || (where 
the matter fields are complex scalar). For Skyrme-Maxwell solitons, toroidal configurations 
were seen only for topological charge two. 

We also looked at the vortex' shape for N = 2 and its dependence on k. The two vortex 
has the shape of a ring for all k, a picture not unfamiliar in planar soliton theories. For 
equivalent coupling, the fields of the two-vortex decay slower than those of the one-vortex. 
This can be understood by looking at formula (|3"8|). Our numerical results show that Voo 
depends strongly on the coupling k but only weakly on the topological charge N such that 
the effective mass k is smaller for the two- vortex and hence its exponential decay slower. 
In Fig. 3 we show the energy density and electric field of the one and two-vortex. The 
coupling here is the lower critical coupling «* r . 
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7 Conclusions 



We have studied classical static vortex solutions in an 0(3)a Chern-Simons system with 
unbroken U(l) gauge symmetry. The vortices have an electric charge which shows a unique 
maximum dependent on the coupling to the gauge field. The magnetic flux in the model 
is effectively quantised in the limit of strong coupling while the angular momentum of the 
vortices is fractional such that they can be considered as classical anyons. 

In the case of two vortices sitting on top of each other, the model has a repulsive and 
two attractive phases, depending on the parameter which couples the gauge and matter 
fields. This has interesting consequences for the interaction of multivortices. In the repulsive 
regime they will presumably try to move away from each other and for a bounded region this 
would lead to a configuration similar to an Abrikosov-lattice with vortices in equidistant 
and fixed positions. Such configurations occur in the description of flux tubes in type- 
II superconductors. In the attractive regime, however, vortices which are not too widely 
separated from each other will be likely to coalesce. In this context it is worth investigating 
whether the vortices of higher winding number show a similar dependence on the coupling, 
in particular whether their critical couplings K l CI and (if they exist) are of the same value 
than they are here. 

The inter-vortex forces at large and medium distances will be dominated by the mat- 
ter fields, because the electromagnetic fields decay faster by a factor of e~ kr . Thus, the 
intercations should be well described asymptotically by the dipole picture developed in [[|. 
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Figure Captions 

Fig. 1 

la) The static energy E (|l6|) in units of 4-7riV as a function of the Chern-Simons coupling 
parameter k for N = 1 (solid line) and N = 2 (dotted line). The plot includes the Bogo- 
mol'nyi bound (dashed line). 

lb) The electric charge Q (|36|) in units of 2ttN as a function of the Chern-Simons coupling 
parameter k for N = 1 (solid diamonds) and N = 2 (triangles up). 

Fig. 2 

2a) The energy density e ( j33|) as a function of r for iV = 1 and k = 0.3 (dotted), 
k = k} ci = 0.632 (solid) and k = 2 (dashed). 

2b) The profile function / as a function of r for N = 1 and «; = «4 r = 0.632 (solid), k = 2 
(dashed), k = 50 (dot-dashed). 

2c) The gauge field a as a function of r for N = 1 and k = = 0.632 (solid), /-c = 2 
(dashed), k = 0.4 (dotted). 

2d) The gauge field v as a function of r for iV = 1 and k = n l cr = 0.632 (solid), k = 2 
(dashed), k = 0.4 (dotted). 

Fig. 3 

3a) The energy density e d33| ) as a function of r for k = re^ r = 0.632, N = 1 (solid) and 
iV = 2 (dashed). 

3b) The electric fields radial component E r as a function of r for «; = k 1 ^ = 0.632, N = 1 
(solid) and iV = 2 (dashed). 
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